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The power sum of P n n components A', , A' 2 , • • • , X„ is defined by 
the relation 

P n = 10 1og in [10 v ' /in + ... + 10 X " /10 J. 

The distributions of such power sums are studied both analytically and by 
Monte Carlo simulation techniques for the case where the components are 
independent, identically distributed, truncated normal random variables. 
Results are given in terms of distributions and moments of P n . The num- 
ber of components varies from 2 to 256, and the standard deviation of the 
component variables before truncation ranges from 1 to 10 dB. The de- 
pendence of the results on the choice of truncation point is also investigated. 

I. INTRODUCTION 

It is common practice in communications engineering to express 
signal and noise powers on a logarithmic scale. As is well known, such 
a scale serves both to narrow the numerical range between large and 
small powers and to simplify some computations by replacing multi- 
plication by addition. The decibel scale is most commonly used. Em- 
ploying this scale, the power level x of a power w is defined by 

w 

x = 10 log 10 - , (1) 

where w is a reference power, and x is expressed in decibels (dB) over 
the reference power w . Note from (1) that w/w = 10 r/1 °. 

In the situation where a number of uncorrelatcd signal sources feed 
into the same load, the power level p„ of a sum of powers u\, . . . , w„ is 
given by 

p„ = 10 1og 10 [10" /,o + ••• + 10- /,0 J, (2) 

2091 



2092 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1967 

where x { is the power level of Wi. Examples of such sums arise in cross- 
talk computations, overload theory for multichannel amplifiers, noise 
calculations on carrier systems and multihop radio systems, and in the 
evaluations of noise distributions on built-up connections between 
telephone subscribers. Here, however, the power levels are in many 
situations random rather than deterministic variables. Thus, in analogy 
with (2) , one is faced with the random variable 

P n = 10 loft. [10 W1 ° + •■• +10 X " /10 J, (3) 

where each X< is a random variable with known distribution. The clas- 
sical power sum problem consists of finding the distribution function 
and the moments of the power sum P„ defined in (3). This problem 
does not, however, possess a simple closed-form mathematical solu- 
tion. As a result, the task of finding approximate solutions has re- 
ceived extensive attention, beginning at least 35 years ago and persist- 
ing till this date. 

Among earlier contributions to the problem, we can distinguish 
those that give specific methods for numerical evaluation of the power 
sum distribution without introducing any other approximations than 
those that are directly related to the numerical technique that is being 
used. 1, 3 - 4 - 5 - 6 Another approach is based on approximating the power 
sum with a normally distributed random variable. 2, 7 This approach, 
due to R. I. Wilkinson, 2 is quite appealing, since it leads to simple 
evaluation formulas. Moreover, it has now been put on a firm mathe- 
matical foundation with the development of a limit theorem by N. A. 
Marlow. In a companion paper, 8 he proves that power sums are asymp- 
totically normally distributed, provided some mild conditions on the 
component variables are satisfied. 

The present paper considers power sums of independent, identically 
distributed, truncated normal random variables, since this is a situa- 
tion of considerable practical importance in transmission engineering 
work. Two approaches are being used. In the first one, asymptotic ex- 
pressions are developed for the mean and variance of P„. The second 
approach is based on Monte Carlo simulation. 9 This method has a 
number of distinct advantages over other numerical methods in that 

(f) it can accept any number of component variables with arbitrarily 

specified distribution functions, 

(u) independence among the component variables is not required, 
(Hi) computation errors do not cumulate as more than two variables 

are added, and 
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(iv) accuracy can be determined through the evaluation of con- 
fidence limits. 

Our main results are numerical estimates of moments of P n and 
selected graphs of its distribution function. A wide range of component 
distributions is covered with n ranging from 2 to 256. Most of the results 
are based on a nominal symmetric truncation of the component variables 
at ±3.5 standard deviations from the mean. In addition, the effect on 
P n of choosing other truncation points is discussed, and some general 
trends are developed. 

II. ANALYTICAL RESULTS 

Consider first the case where the X t are independent, identically 
distributed random variables. Assume that the expectation 



e = E[io x ' /i0 ] 



and the central moments 



t, = E[10 Xl/1 ° - 0]', 

exist and are finite for a sufficiently large range of j. We require — 1 ^ 
; ^ 8 to derive the results for the mean of P n , — 2 22 ; ^ 12 for the 
variance and wider ranges for higher-order moments. 
Rewrite the power sum P n of X x , X a , • • • , X n , as 



P„ = 10 log 10 S n , 



(4) 



where 



Sn = io A " 10 + 



+ 10' 



Now expand (4) in a finite Taylor series about the mean, nd, of S n . 

This gives 



(5) 



Pn 


_ 1 

X 


", , * , S n - nd 1 fS n - nd\ 2 , 
_logM+ nd 2 [ nd ) +••• 

, (-D m + 'AS„ -neY (S n -nd\] 
+ m \ nd J ' lvm \ nO )_ 



where 



X = 



1 



0.23020 



10 log 10 e 

and log stands for log,.. The remainder term in (5) can be expressed in 
integral form as 
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RM -(-DV'jfi^. ,>-i. («) 

or, alternatively, as 

(-1)'" (x 



jtjp) - *=^\-JE—\ , X >-1, (7) 

w m + 1 \1 + 6.17 

where < 8 < 1. 
With B m (x) given by (7) , one obtains 

r J S * ~ nd ) < for m odd, 
\ nd I ' 

so that from (5) we get our first result 

E(P„) ^ LAP n , (8) 

where 

LAP n = 10 log,,, (nd) = | log (nd) (9) 

is the level of average power. 

To derive asymptotic expressions for the moments of P,„ we apply 
the Lemma in Appendix A and (6) to get 

*((^)"W^r/T] " o(«— ) do 

Next, to derive an asymptotic expression for E(P„), we take the ex- 
pected value of both sides of (5) with m = 3. An application of (10) 
then gives 

E{P n ) = LAP n - 5^3 - + 0(l/n 2 ) as n -> qo . (11) 

Here the independence of the component variables has been used to 
express the variance' of S n as nr 2 , and the third central moment of S n 
as nr 3 . The term containing t 3 is of order 1/n 2 . 

To arrive at an asymptotic expression for the variance a*(P n ), we 
use (5) with m = 2 and (11) to get 

P _ E(P) _ 1 S* ~ ne _ A_ fe - ngV 

^" Mn) X n0 2X V nd J 

+ | A ^""/' ) + 0(l/n) as n -> °o . (12) 
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Squaring (12), taking the expected value of both sides, and applying 
(10) to four of the resulting terms gives 



•VJ - frs + <wo 



as n 



(13) 



A similar approach can be used to derive asymptotic expressions for 
higher-order moments. The measures of skewncss and excess, denoted 
by 7i (Pn) and y-2 (P„) , respectively, are defined by 



7,(P„) = 



E(P n - EP„) 3 

*\Pn) 



and 



ftp — pp y 



They are found to satisfy the expressions 
8tT 



Lrf 



-t + 0(l/ri ! ) as n —> co 



and 
y^Pn) = ['I - -^ + -#- - 3_ - + 0(l/n ) as n 



(14) 



oo. (15) 



The asymptotic results given in (11), (13), (14), and (15) are all 
consistent with Marlow's normal limit theorem. 8 The main virtue of 
the asymptotic results above is that they indicate the rate of convergence 
of the four quantities considered. This is of practical interest since 
engineering applications often involve a finite and fairly small number 
of component variables. 

In the particular case where the A', are truncated normal with mean 
dB, standard deviation before truncation of a dB, and symmetric 
truncation at ±c<r dB, the results contained in Appendix B can be 
used to express (11), (13), (14), and (15) in terms of a, c, and n. For 
the mean and the variance we get, respectively, 



i(P n ) = LAP n - 



exp (\V) i- c ((r) -1 

■2\n 



+ 0(l/n 2 ) 



and 



Al\) = f"ZS££VM=± + (l/n 2 ), 

A ll 



(16) 



(17) 
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where the truncation factor U e (a) is defined by 



with 



rr M _ UM 



1 M) ' *(c) - *(-c) 



and 

*(?) = Air f exp(-f/2)d*. 

V27T ''-w 

Derivation of the Wilkinson estimates for the mean and variance of 
the power sum P n is given in Appendix B. This derivation uses the same 
ideas employed by R. I. Wilkinson in 1934. 2 Thus, P n is approximated 
by a normally distributed random variable P nv> . As above, the com- 
ponents are independent, identically distributed truncated normal 
with mean dB, standard deviation before truncation of a dB, and 
truncation at ±c<r dB. From Appendix B we then have 

, (P ..)^^.-5.o g ..[n- CT p< x ' i r- (,r) - i i <»> 



u 



a 



yp ^ 10 , L exp(\V)C/ c (<r) - l l 



(19) 



The first terms in the asymptotic expansion of (18) and (19), re- 
spectively, agree exactly with the results in (16) and (17). This agree- 
ment establishes the important result that expressions (18) and (19) 
are asymptotically correct to the order of n included in (16) and (17). 
Finally, we note that the actual result due to Wilkinson is contained 
in (18) and (19); the case with nontruncated component variables is 
obtained by putting the truncation factor lJ c {a) = 1. 

III. MONTE CARLO RESULTS FOR C = 3.5 

Having established analytical estimates for the mean and variance 
of power sums of truncated normal random variables, let us now turn 
to estimation using the Monte Carlo technique. The power sum problem 
is basically solved by estimating the distribution function of P n . Using 
the Monte Carlo method, one obtains an estimate of this function by 
random sampling. Each sample of the power sum is obtained by selecting 



POWER SUM PROPERTIES 2097 

n independent samples, one from each of the component distributions on 
the dB scale. The corresponding sample value of the power sum is then 
directly computed from (2). For the results presented here, the com- 
ponent samples have been selected via computer generation of so-called 
pseudo-random numbers. These have approximately a uniform distribu- 
tion over the unit interval. Using the inverse error function together with 
nominal truncation at ±3.5<r gave a random variable with truncated 
normal distribution. Because of requirements of computing speed, this 
transformation has been achieved via a table look-up scheme with 
values of the transformation stored in the computer memorj'. 

Table I summarizes Monte Carlo results in terms of estimates of 
the mean /*(P„), the standard deviation a(P„), and the measures of 
skewness and excess yi(PJ and y 2 (P„). Monte Carlo estimates of 
these quantities are denoted by the corresponding latin letters m(P„), 
s{P n ), gi(Pv), and g-j(P n ). The standard deviation and the measures 
of skewness and excess are estimated directly by the corresponding 
characteristics of the sample distribution. The mean is estimated 
through the formula 

m(P n ) = LAP n - (LAP M0 - m MC ). (20) 

The value of LAP n is computed exactly from relation (9), while 
LAPata and m MC are the LAP and the mean, respectively, of the sam- 
ple distribution. The mean fj.(P„) could also be estimated by m MC . 
However, ra(P„) from (20) is preferred over m MC because the Monte 
Carlo results show that it has a smaller sampling variance. 

An indication of the accuracy of the results in Table I is given by 
the number of decimals included. The half-width of the 99 percent 
confidence interval that represents the sampling uncertainty is between 
one and five times the unit in the least significant digit. For the mean, 
the confidence interval width has, however, been computed for m.w C 
instead of for m(P„). The computation of these confidence intervals 
has been based on the asymptotic normality of the corresponding sta- 
tistics. 

Table I shows that the moan of the power sum increases by some- 
what more than 3 dB when the number of component variables is 
doubled for a fixed a. This effect is illustrated in Fig. 1, where the 
mean is plotted as a function of the number of components n. This 
figure shows that the increase in the mean is substantially more than 
3 dB for a doubling of the number of components n in case n is small 
and a is large. On the other hand, Fig. 1 indicates that the slope of the 
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Fig. 1 — Monte Carlo estimates of n(P n ). The components are truncated normal; 
H = 0, truncation at ±3.5(7. 



graph of the mean levels off at approximately 3 dB for each doubling 
of the number of components at all values of o- for n large enough. 

It is illuminating to compare these properties of the mean with the 
properties of LAP,, , According to relation (9), LAP„ increases by 
10 log, „ 2 fti 3 dB for each doubling of the number of components n, 
similar to the increase of the mean noted above. Furthermore, rela- 
tions (8) and (11) imply that LAP n — n(P n ) is nonnegative and ap- 
proaches as n increases toward infinity. The rate of decrease of 
LAP* — n(P n ) is illustrated by the Monte Carlo results plotted in Fig. 2. 

Table I also shows that the standard deviation of the power sum 
decreases as the number of component variables is increased for fixed 
a. This is illustrated in Fig. 3, where Monte Carlo estimates of u(P„) 
are plotted as a function of the number of components n. 

The measures of skewness and excess in Table I can be taken as an 
indication of the deviation from normality of the distribution of the 
power sum. These measures are zero for the normal distribution and 
they have low values for distributions that deviate only slightly from 
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Fig. 2 — Monte Carlo estimates of LAV„-n(P n ). The components are truncated 
normal; n = 0, truncation at ±3. 5a. 
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normality. The table shows that gi and g 2 are very small for rr-values 
up to four over the range of n-values considered. The table also shows 
that g, is, in general, positive. This indicates that the power sum dis- 
tribution is positively skewed. Moreover, g x considered as a function 
of the number of components n has definite maxima around n = 32 
for all sufficiently large values of a. In particular, this means that the 
magnitude of g x decreases as n becomes large enough. This behavior 
is consistent with the asymptotic behavior of the measure of skew- 
ness as expressed by relation (14 1. 

The results of the previous section show that both LAP,, — ^(P,,) 
and a(P„) converge to as n becomes infinite. From these two facts it 
follows that the distribution of P„ — LAP n converges to a distribution 
degenerate at 0. Fig. 4 illustrates this convergence by plots of the 
Monte Carlo estimates of the distribution function of P„ for n = 1,4, 
16, 64, and 256. This convergence is also illustrated in Fig. 5 where the 
1 percent and 99 percent points of the distribution function of P„ are 
plotted in addition to the mean m(P„) and the level of average power 
LAP,,, for a = 10. It is seen that the slope of the 1 percent point with 
a doubling of the number of components can be considerably larger 
than 3 dB, while the 99 percent point changes by somewhat less than 
3 dB whenever the number of components is doubled. LAP n does not 
represent a fixed percentage point on the distribution function as 
n is changing. It is, therefore, seen that the plots in Fig. 5 of some 
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percentage points would actually cross their asymptote LAP n from 
below before approaching it asymptotically from above. This is true 
for all percentage points of the component distribution that lie be- 
tween and LAPi. In other words, the fact that all percentage points 
approach LAP n asymptotically does not imply that the approach is 
monotone. 

IV. COMPARISON BETWEEN ANALYTICAL AND 
MONTE CARLO RESULTS 

At this point it is natural to examine the relative agreement be- 
tween the various analytical approximations and the Monte Carlo 
estimates. Figs. 6 and 7 contain plots of the asymptote (16), the 
Wilkinson approximation (18), and the Monte Carlo estimates of 
LAP n — /j-(P„) for o- = 6 and 10, respectively. Both figures show the 
asymptote as an upper bound for LAP n — ^(P„). The plots also in- 
dicate that the Wilkinson expression gives a better agreement with the 
Monte Carlo results than the asymptote, and they illustrate the de- 
gree of agreement between the Monte Carlo results and the analytical 
expressions for various values of n. Finally, a comparison between the 
two figures shows that the analytical approximations are better fol- 
low values of o- than for high values. Figs. 8 and 9 present similar 
comparisons between Monte Carlo results and analytical approxima- 
tions for a{N„). The figures contain plots of the asymptotic expression 
(17), the Wilkinson expression (19), and the Monte Carlo estimate 
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of a(P„) for <r = 6 and 10, respectively. The figures serve as a basis 
for conjecturing that the asymptote provides an upper bound for <r(P„). 
Furthermore, the figures indicate as above the degree of agreement 
between the analytic approximations and the Monte Carlo results, 
and they show that the analytical approximations are better for low 
than for high values of a. 

V. INFLUENCE OF TAILS 

The results discussed thus far are all based on a truncation of the 
component distributions at ±3.5o\ Truncations at other points can 
easily be studied with the tools used. Thus, Table II summarizes results 
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of Monte Carlo evaluations for symmetric truncations at ±2a, ±2.5<r, 
and ±3o- for a = 1, 6, and 10 dB, respectively, and with the same range 
of ?i-values as considered previously. A study of the table reveals that 
the truncation point can have a considerable influence on the distribu- 
tion of the resulting power sum. To exemplify this, Fig. 10 shows plots 
of the standard deviation of the power sum of 256 components as a func- 
tion of the truncation point c. The plots cover a wider range of c-values 
and rvalues than found in Tables I and II. The extensions are based 
on the Wilkinson approximation. 

The plots in Fig. 10 exhibit the important trend that the influence 
of the truncation point increases with an increase of the component 
standard deviation a. The same conclusion can be drawn from a study 
of the c-dependence of the mean /x(F„) or of the quantity LAP* - 

Table II contains several cases of negative skewness of P„. Hence, 
the earlier observation that P„ is in general positively skewed does not 
apply for c-values below 3.5. 

VI. CONCLUDING REMARKS 

The extension of the results given here to an even larger number of 
components (n > 256) is straightforward, but the computer time 
needed can easily become excessive. The agreement betweent asymp- 
totic expressions and Monte Carlo results for large enough n does, how- 
ever, indicate that the Monte Carlo technique is not necessary for 
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Fig. 10 — Estimates of <r(P2S6)- The components are truncated normal; n = 0, 
truncation at ±ca. 
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power sum evaluations beyond a certain n-value, namely, the one 
where the asymptotic expressions become sufficiently accurate. 

Finally, we note that the problem of evaluating the distribution of 
the power sum of nontruncated normal components has not been 
brought closer to its solution by the results presented here. This prob- 
lem is certainly of mathematical interest even though it represents a 
physically unrealistic situation. Some Monte Carlo studies with larger 
values for the truncation points have indicated that the convergence 
of the power sum to normality is much less rapid in this case, and that 
considerably larger values of the measures of skewness and excess can 
occur than those contained in Table I. 
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APPENDIX A 

Let Xi , X 2 , • ■ • , X„ be independent identically distributed random 
variables. Put 

S n = 10 x ' /in + • • • + 10 x " /10 

and let 6 = E[10 x,/10 ]. In order to prove the asymptotic results in the 
main body of the paper, we need the following. 



Lemma: Suppose 



Q(x) = x' 



,f r r dt T 

Jo 1 + Xt\ ' 



x > -1 



where I, j, m are nonnegative integers. Ij E10 2lx,/1 ° and E10 ,Xl/1 ° are 
bounded, then 



{<<* 



*■■•--« 



nd 



= 0(n~ U2 ) as n -> « . 



Prooj: Let 



r -w - /.' rb dt ' •>-»■ 
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Then Q(x) = rf[I m (x)¥, and it follows from the Cauchy-Schwarz in- 
equality that 

- nd" 



H^r 1 )] \ ^^rW* 



nO 



The asymptotic behavior of the central moment of S n of order 21 is 
found from Cramer. 10 Hence, 



2» 

as n — * <» . 



To complete the proof, it suffices to show that 

■['-P^T - 0(1) -->- 

To show this, we again apply the Cauchy-Schwarz inequality. Thus, 

Jo Jo (1 + zf) 2m + 1 a? + 1 

Hence, 

Consider now the function w(a;) = 1/x', which is convex on (0, °o) for 
j ^ 0. By Jensen's inequality it follows that if «i , • • • , a n ,y y , • • ■ , y„ 
are non-negative real numbers such that «!+•■•+ a n = 1, then 

w(<*i2/i + • • • + a n y n ) ^ ai«(yi) + ■ • ■ + a n u(y H ) . 

In particular, 

(n/S n y = u{SJn) ^ (l/n)[w(10 v,/, °) + • • • + u(10*" /10 )] 

= (l/n)[10- ,Xl/, ° + ••• + lQ- ,Xn/l0 ]. 

Hence, 

E(n/S n y £ E[10-* Xx/1 °]. 

The right-hand side of this inequality is finite by assumption, so the 
proof is complete. 

APPENDIX B 

Derivation of the Wilkinson Results for Truncated Normal 
Components 

As in Appendix A, let X lt X 2 , . . . , X n be independent, identically 
distributed random variables, and assume further that they all have a 
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truncated normal distribution. The density function of Xj is then 

0, x < n — ca, x > m + ca 

(x - tf 



g(.r) = 



I 



•v- 



i I HO V2to CXP 

/I — C(T ^ X ^ H ~\- C(T 

where * stands for the standardized normal distribution function. 

Now let W t be the nonnegative random variable that expresses the 
power corresponding to X it i.e., 

w = 10 x./,o 
The density function of \]\ is 

f0, w < lO "" ,/,0 1 w > I0'" +co,/,n 



}(w) = 



I 



exp 



(log w - Xm) : 



<p(c) - *(-f) \/2tt a\iv '"''V 2XV 

1() („- Cff )/10 ^ w ^ ]() 0, + r,>/.0 

The moments of Wi are therefore, 



'i'TF? =• ( u>*/(to) ^'" = exp (/,-X M + |/c 2 XV 2 )r fl (M, 

•'0 



where 



r«W = 



<l>(c - X<r) - <p(-c - Xa) 



*(c) - <J>(-c) 

accounts for the effect of the truncation. We note that T c {a) — > 1 as 
c — ► 00. 

The mean and variance of ll'i are found to be 



e = E\\\ = exp (X M + |XV)2\(a) 



(21) 



and 



t, - Yar (!!',) = exp (2X/i + XV 2 )r;'(cr)[exp (XV)tf e (<0 - 1], (22) 
where 

n , ? T ,(2a) 

Now let P n be the power sum of X, , A r 2 , • • • , X n and take n = 0. 
Furthermore, let F„ be approximated by a normally distributed random 
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variable P nw . The independence of the X { 's then allows us to establish 
two equations by adding the means and variances of the Wt's to get 
the mean and variance, respectively, of 

S nw = io p "- /l ° = w x + h w n . 

Relations (21) and (22) allow mean and variance of S nw to be ex- 
pressed in terms of mean and variance of P„ w - Hence, we get 

n exp (fAV)!F e (o) = exp \\n(P„„) + £xV(P nte )] 
and 

n exp (\V)7l(er)[exp (XV) U e (o) - 1] 

= exp [2X/x(P„„) + XV(P„„)][exp XV(P n „) - 1]. 
Solving these two equations for /x(P„, r ) and <r 2 (P nw ) we find 



and 



where 



„(P B „) = LAP n - L log [l + ^P (AV)g.(,) - 1 ] 

exp (XV)£/ e («r) - l "| 



<r 2 (P„,) = ^ log 



1 + 



LAP n = i logn + |X<r 2 + i log T e (*). 
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